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We present the redshift-space generalization of the equal-time angular-averaged consistency rela¬ 
tions between {I -|- n)- and n-point polyspectra of the cosmological matter density field. Focusing 
on the case of £ = 1 large-scale mode and n small-scale modes, we use an approximate symme¬ 
try of the gravitational dynamics to derive explicit expressions that hold beyond the perturbative 
regime, including both the large-scale Kaiser effect and the small-scale flngers-of-god effects. We 
explicitly check these relations, both perturbatively, for the lowest-order version that applies to the 
bispectrum, and nonperturbatively, for all orders but for the one-dimensional dynamics. Using a 
large ensemble of Wbody simulations, we find that our squeezed bispectrum relation is valid to 
better than 20% up to IhMpc”^, for both the monopole and quadrupole at z = 0.35, in a ACDM 
cosmology. Additional simulations done for the Einstein-de Sitter background suggest that these 
discrepancies mainly come from the breakdown of the approximate symmetry of the gravitational 
dynamics. For practical applications, we introduce a simple ansatz to estimate the new derivative 
terms in the relation using only observables. Although the relation holds worse after using this 
ansatz, we can still recover it within 20% up to l/iMpc“^, at 2 = 0.35 for the monopole. On larger 
scales, k = 0.2/iMpc~^, it still holds within the statistical accuracy of idealized simulations of volume 
~ 8h“^Gpc^ without shot-noise error. 


I. INTRODUCTION 

Accurate understanding of the nonlinear gravitational 
dynamics is a key for observational projects that mea¬ 
sure the statistical properties of the cosmic structures 
on large scales. The typical scales of interest in these 
projects raMe from the weakly to strongly nonlinear 
regimes [I|, 0. While perturbation theory is expected 
to be applicable as long as the nonlinear corrections are 
sub dominant ii, a fully nonlinear description would 
be helpful to extract cosmological information out of the 
measured statistics over a wider dynamic range. The 
analytical description also becomes more complicated 
when one models higher-order statistics. Although an 
increasing number of analytical techniques to calculate 
the power spectrum or the two-point correlation function 
have been proposed, based for instance on resummations 
of perturbative series expansion or effective approaches 
(sl-fl^. few of them have been applied to the bispectrum 
or even higher orders. 

Consistency relations between different polyspectra are 
then very useful to have an accurate description of the 
higher-order statistics once one has a reliable model for 
the lowest-order one, the power spectrum. Alternatively, 
they can be used to test analytical models, numerical sim¬ 
ulations, or the underlying cosmological scenario (e.g., 
the impact of modified gravity or complex dark energy 
models). Based on the assumption of Gaussian initial 
conditions and gravitational dynamics governed by gen¬ 
eral relativity, these relations hold at the nonperturbative 
level and provide a rare insight into the nonlinear regime 
of gravitational clustering. 


The most generic consistency relations are “kinematic 
consistency relations” that relate the (£ -b n)-density cor¬ 
relation, with £ large-scale wave numbers and n small- 
scale wave numbers, to the n-point small-scale density 
correlation, with I prefactors that involve the linear 
power spectrum at the large-scale wave numbers [13- 
l25l |. These relations, obtained at the leading order over 
the large-scale wave numbers fc', arise from the equiva¬ 
lence principle, which ensures that small-scale structures 
respond to a large-scale perturbation (which at leading 
order corresponds to a constant gravitational force over 
the extent of the small-size object) by a uniform displace¬ 
ment. Therefore, these relations express a kinematic ef¬ 
fect, due to the displacement of small-scale structures 
between different times. This also means that (at this 
order) they vanish for equal-time statistics, as a uniform 
displacement has no impact on the statistical properties 
of the density field observed at a given time. Because 
they derive from the equivalence principle these relations 
are very general and also apply to baryons and galaxies. 
However, in a standard cosmology they provide no infor¬ 
mation at equal times (apart from constraining possible 
deviations from Gaussian initial conditions and general 
relativity). 

To obtain non-vanishing results for equal-time statis¬ 
tics, one must go beyond this kinematic effect. This im¬ 
plies studying the response of small-scale structures to 
non-uniform gravitational forces, which at leading order 
and after averaging over angles correspond to a large- 
scale gravitational curvature. As proposed in [2^ and 
[rTj . this is possible by using an approximate symmetry 
of the gravitational dynamics (associated with the com- 
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mon approximation ft^/f'^ ~ 1, where flm is the matter 
density cosmological parameter and / = d In £>+ /d In a is 
the linear growth-rate), which allows one to absorb the 
change of cosmological parameters (hence of background 
curvature) by a change of variable. These relations again 
connect the {£ + n)-point polyspectra, with £ large-scale 
modes and n small-scale modes to the n-point polyspec¬ 
trum, when an angular averaging operation is taken over 
the £ large-scale modes, which also removes the kine¬ 
matic effect. These consistency relations no longer vanish 
at equal times but they are less general than the previ¬ 
ous relations. Indeed, galaxy formation processes (cool¬ 
ing, star formation,..) introduce new characteristic scales 
that would explicitly break this symmetry. The lowest- 
order relation, which applies to the matter bispectrum, 
has been explicitly tested in using a large ensemble 
of cosmological iV-body simulations, see also [29l - l3ll| for 
related discussions and comparisons with simulations or 
halo models. 


The aim of this paper is to generalize this analysis, 
presented in [2^ and [2^ in real space, to redshift space, 
where actual observations take place. This is only a first 
step towards a comparison with measures from galaxy 
surveys, because we do not consider the important issue 
of galaxy bias in this paper (i.e., to translate our results 
in terms of the galaxy distribution one would need to 
add a model that relates the galaxy and matter density 
fields). However, this remains a useful task as redshift- 
space statistics are well-known to be difficult to model 
because small-scale nonperturbative effects have a non- 
negligible impact up to rather large scales |3, , for 

instance through the fingers-of-god effect [^. Therefore, 
it is even more important than for real-space statistics to 
build tools that hold beyond the perturbative regime. 


This paper is organized as follows. First, in Sec. [II] 
we introduce the statistics of the redshift-space density 
field and its response to the initial conditions. Then, 
in Sec. imi we describe the dynamical equations of the 
system that we consider here and show their symmetry 
that is valid under the approximation — 1. Us¬ 

ing these results, we finally derive the angular-averaged 
consistency relations in Sec. |TVl We focus on the lowest- 
order version of these relations, i.e. the bispectrum, in 
Sec.El where we present our results in terms of the multi¬ 
pole moments of the spectra. We also introduce a simple 
ansatz to estimate new derivative terms in the relation 
from observables, to simplify its form and facilitate the 
connection with practical situations. In Sec. ED the con¬ 
sistency relations are checked both perturbatively and 
nonperturbatively using analytical calculations. We then 
exploit numerical simulations to give a further test of the 
relations in Sec. IVIII We finally summarize our findings 
in Sec. IVIIII 


II. MATTER DENSITY CORRELATIONS 


In this paper, we assume that the nonlinear matter 
density contrast, S{x,t) = [p{x,t) — p\/p, is fully defined 
at any time by the initial linear density contrast 5lo (be., 
decaying modes have had time to vanish), and that the 
latter is Gaussian and fully described by the linear power 
spectrum Pwik), 


((5Lo(fel)^Lo(fc2)) = PLo{ki)SDiki + k2), (1) 

where we denote with a tilde Fourier-space fields. The 
matter density contrast can also be written in terms 
of the particle trajectories, a;(q, t), where q is the La- 
grangian coordinate of the particles, as 


fc 7^ 0 : 5(k, t) 



5{x,t) 


( 2 ) 

(3) 


where we discarded a Dirac term that does not contribute 
for k ^ 0. This expression follows from the conservation 
of matter, pdx = pdq, which yields [1 -I- (5(£c)]d£c = dq. 

Using the Gaussianity of the linear density field 6lo, 
integrations by parts allow us to write the correlation be¬ 
tween £ linear fields and n nonlinear fields in terms of the 
resp onse of the latter to changes of the initial conditions 

mil 


{ho{k[)..5Lo{ki)S{ki,ti)..S{kn,tn)) = 


PLoik[)..PLoiki) 


/ T>^[^(fcl,U)..^(fcn,t„)] \ 

\vho{-k[)..Vho{-k',) / 


(4) 


This exact relation, which only relies on the Gaussianity 
of the initial condition Slo, holds for any nonlinear field 
(5, which is not necessarily identified with the nonlinear 
density contrast. It is also the basis of the consistency 
relations between (£ + n)— and n—point polyspectra, in 
the limit fc' —^ 0, when one can write the right-hand 

side in terms of {6{ki,ti)..S{kn,tn)) multiplied by some 
deterministic prefactors or operators [ITI - I^ . . 

In this paper, we extend the analysis presented in 
for the real-space density field to the redshift-space den¬ 
sity field. Because of the Doppler effect associated with 
the peculiar velocities, the radial position of cosmological 
objects (e.g., galaxies) is not exactly given by their red- 
shift, interpreted as a distance within the uniform back¬ 
ground cosmology. For instance, receding objects appear 
to have a slightly higher redshift than the one associated 
with their actual location, and one is led to introduce the 
redshift-space coordinate s defined as [H, [s^ 


Vr 

s = a; H-Gj., 

a 


(5) 


where is the radial unit vector along the line of sight 
(we use a plane-parallel approximation throughout this 
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article), Vr the line-of-sight component of the peculiar 
velocity v, and a = da/dt the time derivative of the scale- 
factor a{t). Then, the redshift-space density contrast 6^ 
can be written in terms of the Lagrangian coordinate q 
of the particles as [H, [H, [s^ 

fc^O: S^k,t) = I ^ (6) 

where we again discarded a Dirac term that does not 
contribute for k 0. This is the same expression as 
Eq.dS]) for the real-space density contrast S{k), except 
that a;(q, t) in the exponential is replaced by s(q, t), and 
it again follows from the conservation of matter, [1 -1- 
(5®(s)]ds = dq. Then, the redshift-space generalisation 
of Eq.(|3]) reads as 


integral over x' is independent of the position x in the 
small-scale region, at leading order in the ratio of scales, 
and the initial linear density contrast is merely shifted 
by a uniform amount, 

fc' 0 : A6lo = So J dx'W{x')CLo{x'). (12) 

This corresponds to a change of the background mean 
density p, which means that we must obtain the impact 
of a small change of p, hence of cosmological parameters, 
on the small-scale correlation ((5®(si, ti)..(5'*(s„, t„)). 

III. APPROXIMATE SYMMETRY OF THE 
COSMOLOGICAL GRAVITATIONAL DYNAMICS 


(5L0(fc'l)-^L0(W(fcl,il)-5"(fcn,ln)) = 


Pw{k[)..PLo{k'i) 


V5lo{-K)-V5lo{-K) ] 


■(7) 


As in [ 2 ^, [ 2 ^ , we focus on the relations obtained for 
£ = 1 by performing a spherical average over the angles of 
the large-scale wave number k'. This removes the leading 
order contribution, associated with a uniform displace¬ 
ment of small-scale structures by larger-scale modes, that 
vanishes for equal-time statistics (ti = .. = t„) [ITI - I^ . 
One is left with the next-order contribution, which does 
not vanish at equal times [26l - [^ , and is associated with 
the change to the growth of small-scale structures in a 
perturbed mean density background, modulated by the 
larger scale modes. In configuration space, this means 
that we consider angular-averaged quantities of the form 

[Hill 


On scales much smaller than the horizon, where the 
Newtonian approximation is valid, the equations of mo¬ 
tion read as 


|^ + iv[(l + d)^]=0, (13) 

dv 1, , 1 , ^ 

— ^Rv^-{v-\I)v = -V((., 14 

at a a 

= ATiQpa^d. (15) 

Here, we use the single-stream approximation to simplify 
the presentation, but our results remain valid beyond 
shell crossing. Linearizing these equations over {i5, v}, 
one obtains the linear growth rates D±{t), which are the 
independent solutions of 


= Jdx'W{x') {6Lo{x')S‘'{si,ti)..S‘'{sn,tn)), (8) 

which read in Fourier space as 

= {27TfJdk'W{k') {ho{k')5%ki,h)..6%k^,t^)), (9) 

where W{x') [and its Fourier transform W{k')] is a large- 
scale spherical window function. Using Eq. 0, we obtain 


(d«(si,ti)..<5*(s„,U)),„, (10) 


b + 2Hb - AngpD = 0. (16) 

Then, it is convenient to make the change of variables 

[lllilliliS 

ri = lnD+, v = afu, (j) = {dffip, (17) 


with 

d In ab^ 

d In a aZI_|_ ’ 

and the equations of motion read as 


(18) 


and a similar relation for C^, where {■■)eo is the statisti¬ 
cal average with respect to the Gaussian initial conditions 
(5 lo, when the linear density field is modified as 


— -1- V • [(1 -b 6)u] - 0, 

(19) 

Sl{x) ^ SLix) + eoD+{t) J dx'W{x')CLoix,x'), (11) 

i9u , 

^ 2/2 l^u-b(u-V)u- V<p, 

(20) 

where Clo is the linear density contrast real-space cor¬ 
relation function. In the large scale limit for the window 
function W, which corresponds to the limit k' —>■ 0, the 



(21) 
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Here, is the matter density cosmological pa¬ 

rameter as a function of time, which obeys 'i'lrQp = 
As pointed out in [ 2 ^, within the approxi¬ 
mation — 1 (which is used by most perturbative 

approaches [3l), all explicit dependence on cosmology dis¬ 
appears from the equations of motion dm)-dll]). This 
means that the dependence of the density and velocity 
helds on cosmology is fully absorbed by the change of 
variable dni). Then, a change of the background den¬ 
sity, as in Eq.dm, can be absorbed through a change of 
the time-dependent functions f{t)}, which 

enter the change of variables (HZD [ 23 . 

Here, we used the single-stream approximation to sim¬ 
plify the presentation, but the results remain valid be¬ 
yond shell crossing, as the dynamics of particle trajecto¬ 
ries, x{q,t), follow the equation 


d'^x / 3Hni \ dx 

drf y 2p J dr] 


-Vif, 


( 22 ) 


where ip is the rescaled gravitational potential (|2T]) . This 
explicitly shows that they satisfy the same approximate 
symmetry. Therefore, our results are not restricted to the 
perturbative regime and also apply to small nonlinear 
scales governed by shell-crossing effects, as long as the 
approximation flm /~ 1 is sufficiently accurate (but 
this also means that we are restricted to scales dominated 
by gravity). 


Because both frames refer to the same physical system, 
we have r' = r = a'x' = ax, p'(l -|- S') = p{l + S), where 
r = r' is the physical coordinate. Thus, we have the 
relations 

a;'= (1-I-e)a:, <5'= (5 — 3e(l-I-(5), v'= v + eax,{26) 

where we used Ea. (j23l) and only kept terms up to lin¬ 
ear order over e. In particular, we can check that 
if the fields {S',v',(j)'} satisfy the equations of motion 
(HSll-IISl) in the primed frame, the helds {S,v,(f>} sat¬ 
isfy the equations of motion (ll3D - (ll5p in the unprimed 
frame, with the gravitational potential transforming as 
(j)' = (j) — a^{e + 2He)x^/2. This remains valid be¬ 
yond shell crossing: if the trajectories x'{q,t) satisfy the 
equation of motion in the primed frame, the trajectories 
a:(q, t) = (1 — e)a;'(q, t) satisfy the equation of motion in 
the unprimed frame. 

From the dehnition ([5]) and Eq. , we obtain the re¬ 
lation between the redshift-space coordinates, 

S = (1 -1- e)s -t- "rrSr Gr? (^^) 

H 

using a' = {1 — e)a and H' = H — e. Then, using for in¬ 
stance the expressions ([3]) and ([6]), the real-space and the 
redshift-space density contrasts in the actual unprimed 
frame, with the uniform overdensity = 3eo, can be 

written as [ 2 ^ 


IV. ANGULAR-AVERAGED CONSISTENGY 
RELATIONS 

As described in the impact of a small uniform 
change of the matter density background can be obtained 
by considering two universes with nearby background 
densities and scale factors, {p{t),a{t)} and {p'(t), a'(t)}, 
with 

pa^ = pa^=po, a'= a[l-e(t)], p'= p[l-|-3e(t)].(23) 

Here and in the following, we only keep terms up to linear 
order over e. Then, writing the Friedmann equations for 
the two scale factors a{t) and a'(t) and linearizing over 
e, we find that e(t) must sa tisfy the same equation (fTOl) 
as the linear growing mode |38j . Thus, we can write 

eit)=D+{t)eo. (24) 

For our purposes, the universe {p(t),a(t)} is the ac¬ 
tual universe, with the zero-mean initial condition Slq, 
to which is added the uniform density perturbation (1121) . 
To recover zero-mean density fluctuations, we must shift 
the background by the same amount. Thus, this new 
background {p'(t), a'(t)} is given by 

eo = ^ASlo = yJ (25) 

where we used the last relation (P51) . which gives at linear 
order over 5 and e, = <5^ -f 3e. 


5,,{k,t) = ~5[{l-e)k,D+,^], (28) 

and 

= ^®[(1 -e)fc - (29) 

where we disregarded a Dirac factor 5oik) that does not 
contribute for k ^ 0. In Eas. (l28l) - (l2^ . the subscript 
“eo” recalls that we consider the formation of large-scale 
structures in the actual universe to which is added the 
small uniform overdensity ASlo = 3^0 ■ 

The physical meaning of the expression (1^51) directly 
follows from the mapping (1261) and the independence on 
cosmology of the equations of motion (fTOl) - (l?T]) . within 
the approximation ~ 1. It means that in the 

primed universe, with the slightly higher background 
density p' = (1 -I- 3e)p (focusing for instance on the case 
e > 0), comoving distances x' show the small isotropic 
dilatation (1261) [because the higher background density 
yields a higher gravitational force and a smaller scale fac¬ 
tor a'{t)], whence an isotropic contraction of wave num¬ 
bers k', while the linear growth factor D^{t) is also mod¬ 
ified. Moreover, the approximate symmetry discussed in 
Sec. nni implies that all time and cosmological depen¬ 
dence can be absorbed through the time coordinate D^, 
if we work with the rescaled field {(5, u, ip} of Eg. dTTl) . 
This is denoted by the rescaled time coordinate D+jg in 
the right-hand side of Ea. d^ . where the subscript “eo” 
recalls that we must take into account the impact of the 
modihed background onto the linear growth factor. 
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For the redshift-space density contrast (l29t two new 
effects arise, as compared with the real-space density 
contrast (l28)l . First, the mapping s •H’ s' is no longer 
isotropic because of the peculiar velocity component 
along the line of sight, see Ea. (l27l) . which also leads to an 
anisotropic relationship k k'. Second, in addition to 
the time coordinate £>+, the redshift-space density con¬ 
trast involves the new quantity f{t). This follows from 
the definition which can be written in terms of the 
rescaled velocity field u of Eq. (ITTl) as 
Vr 

S = X + — er = X + fUr Gy ■ (30) 

This shows that in addition to the rescaled term Ur, 
which only depends on time and cosmology through the 
time coordinate D+, within the approximate symmetry 
of Sec. mi the line-of-sight component explicitly involves 
a time and cosmology dependent factor f(t), which must 
be taken into account in Ea. d^^ . 

Then, to derive the angular-averaged consistency rela¬ 
tions through Eq. (nni, we simply need to use Eq. (Pll to 
obtain the derivative of the redshift-space density con¬ 
trast with respect to cq, and next to use Eo. d^ . This 
yields 


dS^jk) 

deo 


dD+,, d~S^ dU, d~6^ 

deo dD+ deo df 


dS^ dS^ 


(31) 


where we disregarded a Dirac factor that does not con¬ 
tribute for wave numbers fc ^ 0. 

As found in [2^ |4lj , the derivative of the linear growth 
factor reads as 

= (32) 

£ 0=0 ' 


dP+eo 

deo 


This corresponds to D'^ = D++{13/7)D^eo for the linear 
growing mode in the primed frame, while a' = a — D+aeo 
and H' = H — D+eo- From the definition (fT51) . we obtain 
/' = / + /(13/7D+ -f b+/H)eo, whence 


deo 


£ 0=0 


fD+ 


Therefore, Ea. dHTl) gives 



(33) 


a^"(fc) 

deo 


13 , d5^ 


+ /£>+ I f + / 


, ,. 33' . , 33' 


d5^ 

W 

(34) 


Of course, when we set / to zero, we recover the expres¬ 
sion of the derivative with respect to eo of the real-space 
density contrast S [2^. In configuration space, this reads 
as 


dS^s) 

deo 


D+ 


13 d 
Y dliiD+ 




d d 

+ 3 + f + S ■ — -h fSr -^— 

wS (ys-p 


Sbs) 


(35) 


Next, from Eqs. m. and (1251) . we obtain as for the 
real-space correlations [^ . 


C 7l _ 

vy — 


dx'Wix')CLo{x')J2 


D+i 


i=l 


3 + /i + 


13 d 

Y d\nD+i 





+ fi 




d ' 

F) s 

UOpi _ 


{8bsi,tl)..5{Sn,tn)). 


(36) 


The counter terms of the form — Ijn Sj ensure that all 
expressions are invariant with respect to uniform trans¬ 
lations [by explicitly setting the small-scale region at the 
center of the large-scale perturbation (HU]. They are 
irrelevant for equal-time statistics, ti = .. = tn, where 
factors of the form Si ■ already invari¬ 

ant with respect to uniform translations. 

The comparison with Eq.® gives, after writing the 
correlations in terms of Fourier-space polyspectra. 




{5Lo{k')5\k^M)-Y(kn,tn))'^,^o = PLo{k') 


Itt 

n 


2=1 


1 /, 13 d 


1=1 


n 3n 21 dhiD^i 
n 3 dki 


13 

Y 


+ fi 


fi d 

3 df. 




1=1 


X {^«(fci,ti).. 5 ^(fc„,t„))'. 


(37) 


where ilk' is the unit vector along the direction of k' 
and 6P the Kronecker symbol. The subscript fe' —1 0 re¬ 
calls that this relation only gives the leading-order term 
in the large-scale limit fc' —i 0, whereas the wave num¬ 
bers {fci,.., kn} are fixed and may be within the nonlin¬ 
ear regime. Here we denoted with a prime the reduced 
polyspectra, defined as 

{Sbki).Y{kn)) = {6%ki).Y{k„)y 6D{ki + .. + kr,), (38) 


where we explicitly factor out the Dirac factor associated 
with statistical homogeneity. In particular, this means 
that {S‘^{ki)..S^{kn)y can be written as a function of the 
n — I wave numbers {fci,.., only. 

On large scales we recover the linear theory [3, 
Issjj . with S{k',t') ~ D+(t')Yoik') and S^{k',t') ^ 
D+{t')Yoik'){l -I- /'/i'^), where fi' is the cosine of the 
wave number k' with the line of sight, as in 


k ■ Gr 

M = 


k 


(39) 
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Therefore, Ea. (l37l) also gives 
fd^k' / 


47r \1 + /V2 


6^{ki,ti)..d%K,tr,)} =PL{k\t') 

k' —^0 



- 1 - ./i -(- 

n 3n 21dlnDj^i 

n 

i;(C - 

j=i 

~ r? 3'dkj 

3 j=l 


+ 


y 

1 , h 


+ fi 


fi d 
3 dU 


d 


3 dk 


rj 


In particular, in the linear regime we have the well-known 
expression 


Pl{k;t)=DlPLo{k){l + f^lY, (45) 


where Plo is the linear real-space power spectrum today. 
When / = 0, the relation (l43l) recovers the real-space 
consistency relation, as it should. 




(40) 


When all times are equal, t' = ti = .. = tn = t, this 
simplifies as 


B. Multipole expansion 


dUk' / 6^ik') 


47r \ 1 -I- 


6^{ki)..S^{kn)) =PL{k') 

k'^0 


\ f 13 d fl3 \ / d 

3 21 d\nD+ 7 3 df 


ki d kri d 

^ 3 dki ^ 


i=l 


i=l 


3 dkr 


{S%ki)J^ikn)y. (41) 


V. BISPECTRUM 


The consistency relation (H51) is written for a given 
value of k and /i. In practice, rather than considering the 
redshift-space power spectrum over a grid of one often 
expands the dependence on /i over Legendre polynomi¬ 
als. Thus, we write the nonlinear redshift-space power 
spectrum as 


P%k)=P%k,^l^)=Y,PUk)L2d^k), (46) 

i =0 


A. Relation in {k,iP} space 


The lowest-order equal-time consistency relation ob¬ 
tained from Ea. ddTll corresponds to n = 2, that is, the 
bispectrum built from the correlation between two small- 
scale modes and one large-scale mode. We define the 
bispectrum as in Ea. (l38)) . 


{S%k,)~S^{k,)~S^{ki)) = B%k,,k2,k3) 

x5£)(ki -\- k2 -l- A;3). (42) 


In contrast with the real-space bispectrum, B{ki, k 2 , fca), 
which only depends on the lengths of the three wave 
numbers {ki,k 2 ,k 3 } thanks to statistical isotropy, the 
redshift-space bispectrum also depends on angles because 
the velocity component along the line of sight breaks the 
isotropy. Then, Ea. ddTI) yields 


r drtk' B^{k',k-k'/2,-k-k'/2)k'^o 

J Att 1 -I- 


Pdk') 


' f 13 d 

I -I- - H- 

3 21d\nD+ 




l + /y d 2^,, 2 n d 

3 dink 3^^ ^ ^ ^a/i2 


3df 

p«(fc,y).(43) 


Here we used the symmetries of the redshift-space power 
spectrum to write P®(fc) as a function of k and /i^. In 
Eg. (14311 . the power spectrum is written as a function of 
time through the functions D+ and /, that is, 

P^(fc;t)=P^(fc,y;74+,/). (44) 


where is the Legendre polynomial of order i. Only 

even orders contribute to this expansion because P'* is an 
even function of /r. Substituting into Ea. (l43)) . we obtain 


A£+l 


d/i L2i{iJ.) 


'dOfe. B^{k',k-^,-k-^), 


= PL{k') 

1 d 


f 


l + ^ + 


47r 

13 


d 


1 + 
'13 


PUk) - 


21 d\nD+ 
f d 


7 


+ f 


3df 


3 dink 
{8f+A£-l) 
'(4^-l)(4^-h3) 


{2£-l)2i 

{Ai-3){A£-1) 


P 




3 dInk 


{2i+l){2i+2) 

{A£+3){A£+7>) 


P 


2e+2 


{2£-2){2e-l)2i 


(4£-3)(4£-l) 
{2£+l){2£+2){2£+3) 
{A£+3){A£+b) 


21-2 


21-2 


2i+2 


{2£{2£+l) 

{A£-l){U+3) 

(47) 


For the first two multipoles, 2£ = 0 and 2£ = 2, this 
yields 



'dflk' Bt 


An 1 - 1 - //i'' 


= PL{k') 


-I 


f d Id 




3 df 3 dink 


Poik) 


f d 


3 91nfc 


ipo^(fc) + Ap|(fc) 


13 d 
y ainP 


+ 


(48) 
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and 


which gives 


5 

2 


d/^ L2 (m) 


'dn.' Bt 


f-1 


Air 1 + //x'^ 


= PLik') 


-I 


13 d 


/13 


2151nD+ V 7 


/ 9 15 




3 5/ 3 5 In /c 


/ 5 


3 51n fc 

/ 

3 


11 




21 


^Piik) + ^Piik) 


21 


^2(fc) 


(49) 


C. f-derivative 


1. Relations in space 


In practice, we cannot directly measure the derivative 
with respect to / of the redshift-space power spectrum, 
because the time derivative combines the derivatives with 
respect to -D+ and /. Therefore, the expression (l43l) can 
only be applied to analytical models, where the depen¬ 
dences on and / are explicitly known. To obtain 
an expression that can be applied to numerical or ob¬ 
served power spectra, we must write the derivative with 
respect to / in terms of observed time or space coor¬ 
dinates. Since the redshift-space power spectrum must 
coincide with the real-space power spectrum when either 
/ or vanishes, each factor / (resp. fi^) must appear 
in combination with a power of (resp. /). Here we 
make the ansatz that the dependence on / and pP' only 
appears through the combination which is exact at 
the linear order (|45)l (but at higher orders terms of the 
form fpL^, ..., might appear). This gives 


dHfe, B*(fc', k - k'/% -k - k'/2) ^ ,,,, 

-li - TTTJfi - = 


X 


/ 13 5 1 + /m^ 5 

3 21/5 In a 3 51nA: 



/ + 2/m" 


13a/\ 5 

lap j 3 dp 


P^{k, p;a). (53) 


Using the approximation ilm/P — 1, we might sim¬ 
plify Eq.dSni) by writing / ~ ^[-‘^ + //2 + 3/^/2]. 
However, this introduces an additional source of error, 
and at redshift 2 ; = 0.35, this gives a 15% error on /. We 
checked numerically that this can lead to violations of the 
consistency relations by factors as large as 3 or as small 
as 0.5. Therefore, we keep the expression (1531) in the fol¬ 
lowing. [The impact of the approximation Hm//^ ~ 1 is 
greater on the explicit factor / in Ea. ((53l) than on the 
consistency relation itself, which also relied on this ap¬ 
proximation, because the factor / is evaluated at the ob¬ 
served redshift whereas the consistency relation involves 
the behavior of the growing modes over all previous red- 
shifts, following the growth of density fluctuations, which 
damps the impact of late-time behaviors.] 


2. Multipole expansions 


We can again write the relations (1511) and (1551) in terms 
of the multipole expansion (l46l) . For the first two multi¬ 
poles, Ea. (l5T|l leads to 


P®(fc,t) = P®(fc,/^ ;P+) implies / 


dps 

'W 


= R 


,5P" 


(50) 


This allows us to write Ea. (l43l) as 

f dHfc- B%k', k - k'/2, -k - k'/2) 


47r 


, / 13 

1 + w ' 


1 + /m'2 
5 1-1- ^ p 5 


= PPk!) 


3 2151nP)+ 


3 dink 

P^{kpP-D+).pl) 


In practice, we only measure the dependence of the 
power spectrum with respect to time t, or scale factor 
a(t), and wave number coordinates {fc, p. Then, writing 
d/dt = D+d/dD+-\-fd/df, and using Eg. (1501) . we obtain 


dps 


1 


5 In D 4 


dps 


af 


5 In a af^ dp 


(52) 


d^ f dHfc/ B^, 


k'^O 


l_P J 47r 1 -k fp- 


= PL{k') 


1 5 


3 5 In /c 
65 + 7/ 
210 


P^{k) - ^ 

pm 


/ d 


3 51nfc 
13 + 7/ 
42 


■ / 13 5 

1 + - -I- 

3 2151nP+ 

^Ppk) + Ap|(fc) 






(54) 


and 

IL 


dnL2{p j 

13 5 

'^2151nP^ 

\Poik) 


47r 1 + //r'2 

1 5 


= PL{k') 




3 51n fc 




Pm + Pm 


3 51n/c 
13 + 5/ 


+ 


21 


pm) 


e=3 


126 


(55) 



















































while Ea. (l53l) leads to 


0 


'dnk' Bl, 


fc'-»0 


47r 1 + //j,'- 


1 d 


3 i91n fc 


P^{k) - 


= PL{k') 

f d 


3 9 In /c 

13 / 13a/ \ 

42'^^~42a/2 ) 


' f 13 5 

1 + - H- 

3 21/Ulna 

^o(fc) + ^P2{k) 


/ 13 / 13a/ 

P2{k)+\^+'-^- 


42 6 42a/2 


OO 

xY.P2i{k) 

e^2 


(56) 


and 

u 


dfiL2{fi) 


13 d 


dflfc' 


4 tt 1 + //x '2 

1 d 


= PL{k') 




21/i91na 3 9 In A: 




3 91nfc 


lPoik) + ^P^{k) + ^Pl{k) 


fu 5/ 13a/ 

/65 5/ 65a/^ 

42 + T-W 


'65 65/ 65a/ \ 

^42''"T^~42a/2 j 

OO s 

. (57) 

e=3 ^ 


As compared with Eas. (l4^ and (l4^ . these relations in¬ 
volve all multipoles in the right hand sides, because 
the substitution (l50l) gives rise to factors rather 

than the factor (1 — ^^)d/d^‘^ that appeared in Eq.(|43|). 
In practice, it is not possible to measure or compute all 
multipoles and one must truncate these multipole series 
at some order ^max- This implies an additional approxi¬ 
mation onto these relations O-dSTl). 


VI. EXPLICIT CHECKS 

The angular-averaged consistency relations dSzll-ill]) 
are valid at all orders of perturbation theory and also be¬ 
yond the perturbative regime, including shell-crossing ef¬ 
fects, within the accuracy of the approximation flm//^ — 
1 (and as long as gravity is the dominant process). 

We now provide two explicit checks of the angular- 
averaged consistency relations dSHi-dlll). First, we check 
these relations for the lowest-order case n = 2 , that is, 
for the bispectrum, at lowest order of perturbation the¬ 
ory. Second, we present a fully nonlinear and nonpertur- 
bative check, for arbitrary n—point polyspectra, in the 
one-dimensional case. 


A. Perturbative check 

Here we briefly check the consistency relations for the 
lowest order case, n = 2, given by Ea. (H51) at equal times, 


at lowest order of perturbation theory. At this order, 
the equal-time redshift-space matter density bispectrum 
reads as 0 

H^(fcl,fc2,fc3) = 2DlPLoiki)PLoik2)Ziiki)Zi{k2) 

X Z 2 (fei, fc 2 ) -I- 2 perm. , (58) 


where “2 perm.” stands for two other terms that are 
obtained from permutations over the indices {1,2,3}, and 
the kernels Zi and Z 2 are given by 

Zi(fe) = l + //r2, (59) 


and 

Z2{ki,k2) = 


5 + 3/^^ ^ 1 + fki ^ k2\ 

fei • k2 


7 ' 2 \k2 ' ki J 

kik2 

, ‘2 + dffM^ fki-k2Y , fkfJ' 



7 \ kik2 J ' 2 


X 

^(1 + fk-l) + ^(1 + /Mi) 

, (60) 


where k = ki + k 2 . In the small-fc' limit we obtain 

= 2i^}.Pio(fc')^io(fci)^i(fc')^i(fci) 

X Z2(fc', fci)-I-(fci fc2), (61) 


with fci = k — k'/2 and k 2 = —k~k'/2. Here we used the 
fact that Z 2 {ki, ^ 2 ) vanishes as \ki -l-fc 2 p for \ki +fc 2 | —)■ 
0, whereas PLo{k) ^ fc"® with < 1. [If this is not the 
case, that is, there is very little initial power on large 
scales, we must go back to the consistency relation in the 
form of Eg. (157)) rather than Eg. (HOI) . However, this is 
not necessary in realistic models.] Expanding the various 
terms over k\ as 


Z2{k'M) = 


28 


fe-fc'dPLo 

2k dt + - 

(62) 

k' , „ fc • fc' 

(63) 

fk ■ k'\ ^ 2 

+ 14 M J + + > 


/^'^ k ■ k' k ■ k' k 


+ 


2 kk' 2fc'2 


2 k' 


+ (64) 


substituting into Ea. dHTll . and integrating over the angles 
of fe', we obtain 



1 + /m '2 


PL{k')PL{k){l + 


47 f 

-h - 

21 3 

73/^2 /V' , 4/V^' 
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On the other hand, the right-hand side of Eq. 
at the same order over Pj^ as 


reads 


d^k' Bl 


A-k 1-1- 
1 

2/^2 


= PL{k') 


f 13 d 


l + ^ + 




+ / 


3 9/ 


3 ' 2191nD+ 
1 + d 
3 91n fc 




d 


E>+PLo(fc)(l 


( 66 ) 


Collecting the various terms we can check that we recover 

Eg. (1551). 

Therefore, we have checked the angular-averaged 
redshift-space consistency relation dm for the bispec¬ 
trum, at leading order of perturbation theory, within the 
approximate symmetry flm//^ — 1 discussed in Sec. IIIII 
In this explicit check, the use of this approximate sym¬ 
metry appears at the level of the expression (155)) of the 
bispectrum, which only involves the linear growing mode 
D+ and the factor / as functions of time and cosmology. 
An exact calculation would give factors that show new 
but weak dependencies on time and cosmol^y (and that 
are unity for the Einstein-de Sitter case) [J. These de¬ 
viations from Eq. (l58t are usually neglected [for instance, 
when the cosmological constant is zero, they were shown 
to be well approximated by factors like — 1) that 

are very small over the range of interest [i^]. 

For future use in Sec. I VIIB II below, in terms of the an¬ 
gular monopole and quadrupole. Eg. (1551) gives at lowest 
order of perturbation theory: 

/‘t / + 235/ + 101/> 

,.o.3 ^ Pdk) 35 + 35/ + 21/^ + 5/3dP/.(fc) ] 

' 105 105 dlnfc J ^ ^ 


and 


5 

2 


(420+303/ 

+49/“) Pdk) - § (21 + 18/ + 5/“) (68) 


B. ID nonlinear check 

The explicit check presented in Sec. IVI Al onlv applies 
up to the lowest order of perturbation theory. Because 
the goal of the consistency relations is precisely to go be¬ 
yond low-order perturbation theory, it is useful to obtain 
a fully nonlinear check. This is possible in one dimension, 
where the Zel’dovich solution becomes exact (before 
shell crossing) and all quantities can be explicitly com¬ 
puted. Because of the change of dimensionality, we also 
need to rederive the ID form of the consistency relations. 
We present the details of our computations in App. 
and only give the main steps in this section. 


In the ID case, the redshift-space coordinate ([5]) now 
reads as 


s = X + '^ = X -\- fu, (69) 

a 

where u is the rescaled peculiar velocity defined in 
Eq. (IA4I) , in a fashion similar to Eq. (II3, and the redshift- 
space density contrast ([6l) now writes as 

k^O: 5"(fc,t) = y ^ (70) 

where we again discarded a Dirac term that does not 
contribute for k ^ 0 and q is the Lagrangian coordinate 
of the particles. 

As in the 3D case, to derive the ID consistency rela¬ 
tions we consider two universes with close cosmological 
parameters and expansion rates, a'(t) = [1 — e{t)]a(t). 
Again, from the “ID Friedmann equations” we find that 
e{t) = eoD+(t). Next, a uniform overdensity A^q can be 
absorbed by a change of frame, with eo = ^Slo- Then, 
to obtain the consistency relations, we need the impact 
of the large-scale overdensity A^o on small-scale struc¬ 
tures, which at lowest order is given by the dependence 
of the small-scale density contrast (5®(fc) on eg. As shown 
in App. IA21 this reads as 


dY{k,t) 

deg 




(71) 


As expected, this takes the same form as the 3D result 
([34]) . up to some changes of numerical coefficients. This 
leads to the equal-time redshift-space consistency rela¬ 
tions (see App. IA3p 


y^5®(fci)..5®(fc„)\ =Pi(fcO 

^ J / 

o ^ f) 

+(i+ - Y+f)y^,ki 


1+/+ 


9 




dki 


91nD+ 
(5®(fci)..5®(fc„))'. 

(72) 


Here we no longer need to average over the directions of 
the large-scale wave number k', because at equal times 
the leading-order contribution associated with the uni¬ 
form displacement of small-scale structures by large-scale 
modes vanishes [TtI - I^ . Indeed, because of statistical ho¬ 
mogeneity and isotropy, equal-time polyspectra are in¬ 
variant through uniform translations and cannot probe 
uniform displacements. Therefore, ID equal-time statis¬ 
tics directly probe the next-to-leading order contribu¬ 
tion dUl), which truly measures the impact of large-scale 
modes on the growth of small-scale structures. 

In the ID case, the Zel’dovich approximation is exact 
until shell crossing [2^ and it yields for the redshift- 
space nonlinear density contrast dzni) the expression (see 
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Add.IXH) 

6 ^{k,t)= (73) 

J Stt 

The expression (l7^ is exact at all orders of perturba¬ 
tion theory, but it no longer holds after shell crossing 
(which is a nonperturbative effect). On the other hand, 
we can define a ID toy model by setting particle trajec¬ 
tories as equal to Ea. (IA17l) . This system is no longer 
identified with a ID gravitational system, and it only co¬ 
incides with the latter in the perturbative regime, but it 
remains well defined and given by Eas. (IA17p and (l73l) in 
the nonperturbative shell-crossing regime. 

Then, using the expression (1751) we can explicitly check 
the ID consistency relations d72|) . We present in App. IASI 
two different checks. Eirst, in App. lA 5 al we check 
Eg. (1711) by explicitly computing the impact on the non¬ 
linear density contrast (1751) of a small change A6lo to 
the initial conditions. Second, in App. lA 5 bl we directly 
check the consistency relations (1751) by explicitly com¬ 
puting the correlations (^L(fc')(5®(fci)..(5®(fc„))).,^Q and 
(5®(fci)..5'*(fc„))' and verifying that they satisfy Eg. (1751) . 

These two different checks allow us to check both the 
reasoning that leads to the consistency relations, through 
the intermediate result m, and the final expression of 
these relations. They also explicitly show that they are 
not restricted to the perturbative regime. In particular, 
they extend beyond shell crossing, as seen from the toy 
model defined by the explicit expression (175)) (i.e., where 
one defines the system by the Zel’dovich dynamics, even 
beyond shell crossing, without further reference to grav- 

ity)- 

As for the real-space consistency relations , it hap¬ 
pens that in this ID model (1751) the ID consistency rela¬ 
tions dHI) are actually exact, that is, they do not rely on 
the approximation k ~ kq; where k defined in Ea. dASI) 
plays the role of the 3D factor encountered in 

Eos. (fT51) - (15T1) . This is because the redshift-space density 
contrast (1731) truly only depends on cosmology and time 
through the two factors and /, even at nonlinear or¬ 
der. In contrast, in the 3D gravitational case, beyond 
linear order new functions of cosmology and time appear 
(for cosmologies that depart from the Einstein-de Sitter 
case) and they can only be reduced to powers of D+ and 
/ within the approximation — 1- On the other 

hand, if we consider the actual ID gravitational dynam¬ 
ics even beyond shell crossing, where it deviates from the 
expression dza, then the ID consistency relations (1751) 
are only approximate in the nonperturbative regime, as 
they rely on the approximation k ~ kq, while remaining 
exact at all perturbative orders. 

Unfortunately, it is not easy to build 3D analytical 
models that can be explicitly solved and suit our pur¬ 
poses. The 3D Zel’dovich approximation again provides 
a simple model for the formation of large-scale structures 
and the cosmic web. However, it cannot suit our purposes 
because it does not apply to the dynamics of the 3D back¬ 
ground universe itself. Indeed, as can be seen from their 


derivation in Sec. IIVI the consistency relations precisely 
derive from the fact that a large-scale almost uniform 
density perturbation can be seen as a local change of the 
cosmological parameters (i.e., the background density). 
This is also apparent through the fact that the devia¬ 
tion e{t) between the two nearby universes (1551) obeys the 
same evolution equation (1161) as the linear growing mode 
of local density perturbations. This is no longer possi¬ 
ble for the 3D Zel’dovich approximation, which is not 
an exact solution and cannot be extended to the Hubble 
flow itself. In contrast, in the ID universe the Zel’dovich 
approximation is actually exact (before shell crossing) 
and it applies both at the level of the background and 
of the density perturbations. An alternative dynamics, 
which is exact at the background level and provides an¬ 
alytical results on small nonlinear scales, is the spherical 
collapse model. However, this yields a very different den¬ 
sity field than the actual one, as there is a single central 
density fluctuation that breaks statistical homogeneity 
and density correlations are no longer invariant through 
translations. Therefore, although it should be possible 
to obtain some consistency relations for this model, they 
would have a rather different form and this ID spheri¬ 
cal model would be even farther from the actual universe 
than the ID statistically homogeneous model studied in 
this section. 


VII. SIMULATIONS 

The angular-averaged consistency relations (1571) - (15T1) 
are valid at all orders of perturbation theory and also be¬ 
yond the perturbative regime, including shell-crossing ef¬ 
fects, within the accuracy of the approximation — 

1 (and as long as gravity is the dominant process). We 
have explicitly confirmed them either perturbatively or 
nonperturbatively, but the latter is limited to the one¬ 
dimensional case. 

It would thus be of great importance to further check 
these relations in three dimensions nonperturbatively. 
We here exploit a series of A-body simulations for this 
purpose. As can be seen in the following, they are also 
useful to understand the possible breakdown of the rela¬ 
tions and test the validity of the ansatz employed in the 
measurement in practical situations. We first summarize 
how we can evaluate the derivative terms in the consis¬ 
tency relations. We then present the numerical results 
for the bispectrum together with a brief description of 
the simulations themselves. 


A. Derivatives from numerical simulations 

The consistency relation (1551) involves derivatives with 
respect to D+ and /. They can be obtained at once 
within the framework of an explicit analytic model for 
the matter density polyspectra. However, in this paper 
we do not use these relations to check a specihc analytical 
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model. Instead, we wish to use numerical simulations to 
test these relations (which are only approximate because 
of the approximation — !)■ Nevertheless, we can 

also measure separately the derivatives with respect to 
Z1+ and / from the simulations. 

The redshift-space coordinate s can be written in terms 
of the comoving coordinate x and peculiar velocity v as 
in Eg. (I, SOI) . As explained in Sec. IIIII within the approxi¬ 
mation — 1 that is used to derive the consistency 

relations, all time dependence can be absorbed in the 
linear growing mode D^(t) with the change of variables 
(HU. This means that the fields {J, u, </?} are only func¬ 
tions of time through D+, as well as the displacement 
field (q, t) = a; — q, where q is the Lagrangian coordi¬ 
nate of the particles. Thus, for a given realization defined 
by the linear density field ^Lo(q) (normalized today or at 
the initial time of the simulation), the redshift-space co¬ 
ordinate s depends on the functions D+{t) and f{t) as 

s{q,t) = x{q,D+) + fur{q,D+)er. (74) 

Therefore, a small change A/ of the factor / corresponds 
to a change of the redshift-space coordinate s(q) of the 
particles given by: 

/—>■/ + A/ : s —s -I- A/ UrBr = s + A/^ e^. (75) 

On the other hand, from the equations of motion (fT^ - 
(EH, a change A In D+ of the linear growing mode leads 
to a change of the particle velocities and coordinates 


In In £>+ -f- A In £)+ : a; — >• a: + A In £)+ u, 

/3f2, 


u — >• u — A In £) 

whence, 

s — >■ s -I- A In £)+ u — / 
= s 4- A In D+ ^ - 

[af 


+ \2p 

I 2/2 

3n„ 


-1 u + Vy^ , (76) 


1 ) Ur „ 1 Gt- 

or 


.Vr 1 dc/)] 


(77) 


Thus, to obtain the partial derivative of the power 
spectrum with respect to / or lnll+, we modify the par¬ 
ticle redshift-space coordinates by Eqs. (ESI or ([77|. for a 
small value of Af or A In!?_)_, and we compute the as¬ 
sociated power spectrum. Taking the difference from the 
initial power spectrum and dividing by Af or A In 
gives a numerical estimate of dP^/df or dP^/dliiD^. 


B. Numerical results 

We are now in a position to present the consistency 
relations measured from simulations. Before that, let 
us briefly describe the simulations used here. They are 


the ones performed in Q. Employing 1024^ dark mat¬ 
ter particles in a periodic cube of (2048 h“^Mpc)^, the 
gravitational d yna mics was solved by a public simulation 
code Gadget 2 starting from an initial condition set 

at 2 = 15 by solving second-order Lagrangian perturba¬ 
tion theory |4^ |4^. The cosmological model used 
was a flat-ACDM model consistent with the hve-year 
observation of the WMAP satellite [d^: f2„i = 0.279, 
fib = 0.165flm, h = 0.701, As = 2.49 x 10“® and 
rig = 0.96 at fco = 0.002Mpc~^. This whole process was 
repeated 60 times with the initial random phases varied 
to have a large ensemble of random realizations. 

The consistency relations have already been examined 
and presented in real space in [^ . There it was found 
that the relation was recovered within the numerical ac¬ 
curacy at z = 1, while discrepancy of several percent 
level was found a,t z = 0.35. It was further discussed 
that this is presumably due to the breakdown of the ap¬ 
proximation Om//® ~ 1; we could indeed confirm that 
the relations better hold in supplementary simulations 
done in the EdS background, but with exactly the same 
initial perturbations. We focus here on the lower red- 
shift, z = 0.35, at which the consistency relations are the 
most nontrivial. 


1. Full consistency relations 

We first consider the redshift-space consistency rela¬ 
tions in their full form (H51) - (H^ . with both derivative 
operators d/dD^ and d/df. 

Having already presented the methods we employ to 
measure the derivative terms in the previous subsection, 
the post-processing for the simulation outputs is exactly 
the same as in except that we now consider the par¬ 
ticle positions in redshift space. The matter density field 
is constructed with the Cloud-in-Cells (CIC) interpola¬ 
tion on 1024® mesh cells and subsequent computations 
are based on the fast Eourier transform. The change in 
the particle coordinates corresponding to a slight change 
in In £)+ is also computed based on the calculation on the 
same mesh cells for d(j)/dr and then interpolated to the 
positions of particles using the CIC kernel (see Eq. [TZl). 

The monopole and the quadrupole moments of the re¬ 
lation for the bispectra, Eas. (H51) and (03), are respec¬ 
tively shown in the left and the right panel of Fig. [T] 
In each panel, we fix the value of the hard wave mode 
k, and plot the ratio of the two sides as a function of 
the soft mode k'. The error bars are estimated based on 
the scatter among the 60 independent realizations. They 
thus correspond to the error level expected for an ideal 
survey with a volume of ^ 8h“®Gpc® when we can ig¬ 
nore the shot noise contamination. Overall, the ratio is 
close to unity for both monopole and quadrupole. From 
this hgure, we basically confirm the relations at the non- 
perturbative level in the three dimensional dynamics. 

The dashed lines in Fig. [T]show the ratio of the mea¬ 
sured bispectrum to its tree-order predictions (1671) and 
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FIG. 1; Consistency-relation ratio for the redshift-space bis- 
pectrum from A'^-body simulations. The dashed lines show the 
ratio of the measured bispectrum to its tree-order predictions 
(|67)l and ((Ml)- 


(IB51) . For the monopole, this lowest-order perturbative 
prediction fares reasonably well as it only underestimates 
the nonlinear results by 30%, on these scales. How¬ 
ever, it is already less accurate than our result (H51) . 
which takes into account higher-order and nonperturba- 
tive nonlinear corrections (at the price of the approxima¬ 
tion ~ 1). For the quadrupole, the lowest-order 

perturbative prediction does not appear in the panels 
at fc > 0.4/iMpc“^ because in these cases it is out of 
range and actually gives the wrong sign. This change 
of sign is likely due to the fingers-of-god effect, which is 
not captured by perturbation theory. Indeed, it is well 
known that higher-order multipoles are increasingly sen¬ 
sitive to small-scale nonlinear contributions, as finger-of- 
god effects immart a strong angular dependence to the 
bispectrum |32| . In contrast, our result (1491) remains 
consistent with the simulation data within 20%. This 
shows that we test the consistency relations in a nontriv¬ 
ial regime, beyond the reach of standard perturbation 
theory. Thus, the trade-off between the error introduced 
by the approximate symmetry of Sec. IIIII and the advan¬ 
tage of taking into account all nonlinear contributions, at 
both perturbative and nonperturbative levels, is benefi¬ 
cial. This is particularly true for complex statistics such 
as the redshift-space quadrupole that are very sensitive 
to small-scale highly-nonlinear effects, which are difficult 
to include in analytical modelings. 

However, when we look into each panel more closely, we 
can find that the data points are slightly off from unity. 
For the monopole moment, the ratio tends to be larger 


than unity at A: 0.4hMpc“^. On the other hand, unity 
is within the statistical error level for the quadrupole mo¬ 
ment, though the central values are larger (smaller) than 
unity on k ^ 0.4/iMpc“^ (fc ^ 0.8 fcMpc”^). In most of 
the cases, the deviation from unity is at most 20%, and 
this is meaningful only when we measure the ratio very 
precisely; an ideal survey with a volume of ^ 8 /i“^Gpc^ 
can detect the deviation from unity only for the monopole 
moment on small scales. 

These deviations are somewhat greater than those 
found in in real-space, which only reached 7% at 
fc = IfcMpc”^. This is not surprising, because it is well 
known that redshift-space statistics are more sensitive to 
small nonlinear scales, for instance through the fingers 
of god effect, and low-order perturbation theory has a 
smaller range of validity. Then, we can expect a greater 
violation of the redshift-space consistency relations be¬ 
cause the breakdown of the approximation ~ 1 

has a stronger impact on higher perturbative orders. In¬ 
deed, absorbing the time and cosmological dependence 
by D^{t) and f{t) is exact at linear order whereas higher 
orders involve new functions (t) that are not exactly 
equal to and the discrepancies may cumulate 

in the nonlinear regime. 



FIG. 2: Same as Fig. [T] but for the simulations performed in 
EdS background expansion. 

We then work on the supplemental simulations done 
in the EdS background, to understand the cause of this 
small discrepancy, just as in our previous real-space paper 
[ 2 ^. Note that our consistency relations in an Einstein- 
de Sitter cosmology also involve the approximate sym¬ 
metry described in Sec. uni even though = 1 

in the EdS background. Indeed, what matters is not 
that ^ra/P be unity in the reference cosmology, but that 
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^m//^ remain (approximately) constant as we vary the 
background curvature around the reference cosmology. 
Nevertheless, the comparison between EdS and ACDM 
results provides a simple estimate of the impact of our 
approximation, because the difference between these two 
cosmologies arises from the change of reference point 
along the ilmlP curve. 

The results from four realizations of such simulations 
are shown in Fig. [2] Although the scatter of the data 
points are larger than in Fig. [l] the systematic depar¬ 
ture from unity in the previous figure is clearly reduced. 
We thus conclude that the small violation of the consis¬ 
tency relations for the bispectrum can be explained by 
the breakdown of the approximation — 1 (more 

precisely, of constant for nearby background cur¬ 

vatures), in agreement with the discussions above. 

2. ffi^ ansatz and reduction to d/dD+ operator 


infinite summation over all Legendre multipoles of the 
redshift-space power spectrum. Here, we truncate these 
series at order (crosses) or PI (pluses). 

The difference between the two symbols is negligi¬ 
ble at fc = 0.2 and 0.4/iMpc“^ for the monopole and 
fc = 0.2 /iMpc”^ for the quadrupole moment, where the 
ratio itself is roughly consistent with unity. As we move 
to smaller scales, the two symbols become more distinct. 
In those cases, adding the higher-order term (i.e., P4) 
does not help to restore the relations, suggesting that the 
ansatz is not a good approximation at the corresponding 
scales. The plot suggests that the quadrupole moment 
is more sensitive to the higher-order term and thus the 
ansatz works less accurately than for the monopole mo¬ 
ment. This is naturally expected since the quadrupole 
moment is impacted more strongly by higher-order cor¬ 
rections [see e.g. (48l |. where we can see how much higher- 
order perturbative corrections affect the first two mo¬ 
ments. These corrections have terms where m 

and n can be different.]. 



FIG. 3: Same as Fig. [T] but for Eqs. (1541) and (1551) after apply¬ 
ing the ansatz to rewrite the /-derivative terms. Blue crosses 
and red plus symbols respectively show the ratio truncated at 
Pi and Pi. 

Now, we come back to the original ACDM simulations 
and apply the ansatz that all the / and ^ dependences ap¬ 
pear through the combination //i^. This allows us to re¬ 
place the /-derivatives, djdf, by /r^-derivatives, 3/9/i^, 
as in Eg. (1501) . This gives the approximated consistency 
relations for the bispectrum, Eqs. (IMl) and (l55l) . respec¬ 
tively for the monopole and quadrupole moment, which 
we plot in Fig. [3] In contrast with the exact form of 
the consistency relations, given by Eas. lHHl) and (H3|) and 
displayed in Fig. [U the right-hand side now involves an 


3. fjP ansatz and further reduetion to d/da operator 
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FIG. 4: Same as Fig. [S] but for Eqs. (1561) and (1571) after 
applying the ansatz once again to convert the D+-derivative 
terms into the a-derivatives. 

The situation is basically the same after we further 
apply the ansatz to replace the derivative with respect 
to D+ by a derivative with respect to time or the scale 
factor, as in Ea. (l5^ . As compared with the form dis¬ 
played in Fig. [3l this involves an additional approxi¬ 
mation, which relies on the same //i^-ansatz, because 
the full time derivative, or scale-factor derivative, d/da, 
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combines both theoretical derivatives d/dD+ and djdf. 
Therefore, to replace the operator d/dDj^ by d/da we 
must once again use the //i^-ansatz to remove the new 
d/df terms generated by the change of variable from 
to a. 

Figure |4] shows the results of Eqs. (l56l) and (l57l) . The 
plus symbols are now out of the plotted range for the 
quadrupole moment on small scales. The relation for the 
monopole moment is more robust against this approxima¬ 
tion ansatz on large scales, especially at fc = 0.2 hMpc”^, 
and we can safely apply the consistency relation here in 
the simplified form (1551) . Except for this case, the ratio is 
affected significantly by the ansatz and the order at which 
we truncate the infinite summation on the right-hand 
side. Nevertheless, we note that by truncating at order 
we obtain a good agreement, better than 20% up to 
k = l/iMpc“^, for the monopole. For the quadrupole, 
the deviation can reach up to 40%. 

Therefore, even with the current ansatz, we can still 
examine how the ratio behaves in the observations and 
compare it with the simulation results. Since the data 
points obtained with the truncation at order (i.e., 
cross symbols) are less noisy and moreover stay around 
unity after applying the ansatz, the easiest check of the 
true gravitational dynamics is to apply the same ansatz 
and truncate the moments at this order. We would need a 
more involved ansatz for the estimation of the derivative 
terms from observations to extend the applicable range of 
these consistency relations, and we leave this to a future 
study. 


VIII. SUMMARY 

In this paper, we have generalized the equal-time 
angular-averaged consistency relations for the cosmic 
density field originally developed in real space by [2^ 
to redshift space, in which the actual observations are 
taking place. These relations express the squeezed limit 
of (n-l-l)-point correlation functions or polyspectra, with 
n small-scale modes (that can be in the nonlinear regime) 
and one large-scale mode (in the linear regime at a much 
larger scale than all other n wave numbers), in terms of 
the n-point correlation of the small-scale modes. These 
relations can be generalized to (n + £) correlations, with £ 
large-scale modes, as in [2^, but we focused here on the 
case of one large-scale mode. The explicit forms that we 
have obtained rely on an approximate symmetry of the 
dynamics, — 1- However, within this approxima¬ 

tion they are valid at a fully nonlinear level. Thus, they 
hold at all orders of perturbation theory and also in the 
nonperturbative regime, beyond shell crossing. In partic¬ 
ular, they include both the large-scale Kaiser effect 
associated with the infall of matter within large-scale 
gravitational wells, and the fingers-of-god effect [^, as¬ 
sociated with the virial motions inside collapsed halos. 

We have found that, because the mapping from real 
to redshift space involves the velocity component along 


the radial direction, the form of these consistency re¬ 
lations is slightly more complex than in real space, as 
it involves two types of time derivatives. The first is a 
derivative with respect to the linear growing mode (t ), 
which also appeared in the real-space case. The second 
is a derivative with respect to the linear growth rate, 
f{t) = dlnZ)+/dlna. This differential operator, d/df, 
did not appear in the real-space case and it arises from 
the scaling of the peculiar velocity field (i.e., through the 
change of variable from to u, where u is the rescaled 
velocity field that makes use of the approximate symme¬ 
try of the dynamics). This feature makes it more difficult 
to use these relations for observations, because at best we 
can only measure one time derivative, d/dt, which com¬ 
bines both d/dD+ and d/df, and we cannot measure 
these two derivatives separately. However, these relations 
can still be used to check analytical models or numerical 
simulations, where we can explicitly compute these two 
derivatives. 

Next, we have tested these consistency relations both 
analytically and numerically. First, at leading order of 
perturbation theory, we have checked the lowest-order 
consistency relation, which expresses the squeezed limit 
of the bispectrum in terms of the nonlinear power spec¬ 
trum of the small-scale modes. Second, in a fully non¬ 
linear and nonperturbative analysis, we have checked 
all these consistency relations at all orders, in the sim¬ 
pler one-dimensional case, where we can use the exact 
Zel’dovich solution of the dynamics. 

We have also tested the lowest-order consistency rela¬ 
tions, relating the nonlinear bispectrum and power spec¬ 
trum, with numerical simulations. We find a reasonably 
good agreement at z = 0.35. Projecting the angular de¬ 
pendence of the redshift-space polyspectra onto Legendre 
polynomials, we find a good agreement for the monopole 
up to fc ;< 0.4ft,Mpc“^ and we detect a small deviation of 
at most 20% for k < l/iMpc“^. For the quadrupole, we 
do not detect significant deviations (but the statistical 
error bars are slightly larger). In the case of an Ein¬ 
stein de Sitter cosmology, we find that these deviations 
are greatly reduced and our numerical data agree with 
theoretical predictions. Therefore, the small deviations 
found in the ACDM cosmology can be explained by the 
finite accuracy of the approximation — 1- 

The typical magnitude of these deviations is larger and 
extends over a wider wave number range than for the real- 
space consistency relations [2^. This is consistent with 
the observation that the nonlinearity in the cosmic veloc¬ 
ity field is more sensitive to the local nonlinear structure 
on small scales, as small-scale effects can easily propa¬ 
gate to larger scales through the nonlinear mapping from 
real to redshift space. Indeed, it is well known that the 
perturbation-theory prediction of the matter power spec¬ 
trum is more difficult in redshift space 0 . Then, because 
nonlinear effects are likely to amplify the breakdown of 
the approximation ~ 1, violations of the consis¬ 

tency relations due to breakdown of this approximate 
symmetry are indeed expected to be greater in redshift 
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space. 

On the other hand, we find that our results for the bis¬ 
pectrum provide a significant improvement over lowest- 
order perturbation theory, especially for the quadrupole 
where the perturbative prediction even gives the wrong 
sign for k > 0.4/iMpc“^. This is a signature of the strong 
impact of small-scale nonlinearities onto redshift-space 
statistics, which are usually difficult to model analyti¬ 
cally. This shows that we test the consistency relations in 
a nontrivial regime, beyond low-order perturbation the¬ 
ory. It also shows the interest of these nonlinear relations, 
as the inaccuracy introduced by the approximate symme¬ 
try ~ 1 is more than compensated by the account 

of higher-order and nonperturbative nonlinear contribu¬ 
tions. This can be even more beneficial for statistics 
such as the redshift-space quadrupole that are sensitive 
to highly nonlinear effects that are difficult to model. 

To make the connection with observations, or to sim¬ 
plify the form of these consistency relations, we also 
tested a simple ansatz that allows to remove the new op¬ 
erator djdf. This relies on the approximation that / and 
/i^ only enter the redshift-space power spectrum through 
the combination (this is exact at linear order, in the 
Kaiser effect). A first step allows us to remove the oper¬ 
ator d/df, which only leaves the operators d/dD+ and 
d/dk in multipole space. The drawback is that the right- 
hand side of each consistency relation now involves an 
infinite series over multipoles of all orders. We find that 
this approximation gives rise to an additional source of 
discrepancy between the numerical data and the analytic 
predictions, especially for the quadrupole. Moreover, the 
result depends on the order at which we truncate the 
multipole series in the right-hand side. It turns out that 
better results are obtained when we truncate at the low¬ 
est order P|. This suggests that the //r^-ansatz does not 
faithfully describe higher perturbative or nonperturba¬ 
tive orders. 

In a second step, we use once more the //r^-ansatz to 
replace the operator d/dD+ by the full time-derivative, 
or scale-factor derivative d/da. As could be expected, 
we find that this further increase the deviation from the 
numerical simulation data and the dependence on the 
truncation order, especially for the quadrupole. 

Nevertheless, we find that using the /^^-ansatz and 
truncating at order PJ* we obtain an agreement that is 
better than 20% up to fc = l/iMpc“^ for the monopole. 
For the quadrupole the deviation can reach 40%. Al¬ 
though these limitations make the accessible range of 
these consistency relations rather narrow, we can still 
make use of them to predict the higher-order polyspec¬ 
tra. There is substantial recent progress on the redshift- 
space clustering, but the calculations are mostly limited 
to the power spectrum (e.g., [^ Idsl - ls^ l. Using the rela¬ 
tions developed here, one can compute, for instance, the 
angular-averaged bispectrum in redshift space by substi¬ 
tuting these formulae for the power spectrum. Since the 
relations approximately hold down to very small scales, 
albeit not perfectly, they can be useful in estimating the 


covariance matrices of the redshift-space observables [we 
need the trispectrum to compute the matrix for the power 
spectrum]. Indeed, the accuracy required for the covari¬ 
ance matrices might not be as demanding as that for the 
spectra themselves. A study along this line is undergoing 
now, and we wish to present the results elsewhere in near 
future. 

A more complex issue is the problem of biasing, when 
we wish to connect measures from galaxy surveys with 
theoretical predictions. In principle, the approximate 
symmetry flm/— 1 that we used to obtain explicit 
expressions no longer applies once we take into account 
galaxy formation physics. Indeed, baryonic processes 
(cooling, star formation, ....) involve new characteristic 
scales that explicitly break the symmetry of the dynam¬ 
ics. Then, a priori it is no longer possible to absorb the 
time-dependence of the dynamics by a simple rescaling 
that only involves the linear growing mode. Therefore, 
the relations we have obtained are not guaranteed to ap¬ 
ply to the galaxy density field itself, by making the naive 
replacement S ^ Sg. One should rather use these rela¬ 
tions as constraints on the matter density field, and given 
a supplementary model that relates the galaxy field to the 
dark matter density field, derives the consequences onto 
the galaxy density field. Of course, this would depend on 
the model that is used to describe galaxy formation and 
introduce an additional approximation. We leave such a 
study for future works. 
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Appendix A: ID example 

As for the real-space consistency relations [2^ , it is in¬ 
teresting to check the redshift-space consistency relations 
(l37l) - (l4T]) obtained in this paper by using a simple one¬ 
dimensional example that can be exactly solved. This is 
again provided by the Zel’dovich dynamics (4^ , which is 
exact in ID (before shell crossing). 

1. ID equations of motion 

The ID version of Eos. (fT^ - (ITS]) reads as 
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dv 1 dv 1 d(j) 

+ Hv+ -V— = 


(A2) 


^ = A^g{t)-paH. (A3) 

Here we generalized the ID gravitational dynamics to the 
case of a time-dependent Newton’s constant Q{t). This 
allows us to obtain ever-expanding cosmologies, similar 
to the 3D Einstein-de Sitter cosmology, for power-law 
cases Q{t) oc with — 2 < a < — 1 [and a{t) oc 
p{t) oc l/a(<) oc 

Linearizing these equations, we obtain the evolution 
equation of the linear modes of the density contrast. It 
takes the same form as the usual 3D equation (fTO . D + 
2H{t)D — 4:'KQ{t)p(t)D = 0, but with a time-dependent 
Newton’s constant and the ID scale factor a{t). 

In a fashion similar to the change of variables (ED , we 
make the change of variables 


r] = \nD+, v = afu, (j) = {af)'^(p, with / 


aii)+ 

dD+’ 


(A4) 


and we obtain the rescaled equations of motion 

| + A[(1 + «„|^0. (A5) 


du , ,, _ du 

^11“+“& 


d(p 


(A6) 


^ = K{t)6, (A7) 

where we introduced the factor K{t) defined by 

n{t)=4ngit)p{t)^^- (AS) 

Thus, K{t) plays the role of the ratio 3Dni/(2/^) encoun¬ 
tered in the 3D case in Eas. (fT^ - (|21l) . Then, the 3D 
approximation — 1 used in the main text corre¬ 

sponds in our ID toy model to the approximation k ~ kq. 
That is, we neglect the dependence of k on the cosmo¬ 
logical parameters and time, and the dependence on the 
background is fully contained in the change of variables 
(IA4I) . [The generalization to the case of a time-dependent 
Newton’s constant is not important at a formal level, be¬ 
cause it does not modify the form of the equations of 
motion. However, it is necessary for this approximate 
symmetry to make practical sense, so that we can find 
a regime where k is approximately constant. This cor¬ 
responds to cosmologies close to the Einstein-de Sitter¬ 
like expansion a{t) oc in the case g{t) oc with 

-2 < a < -1.] 

The fluid equations (IA5D-(IA7D only apply to the single¬ 
stream regime, but we can again go beyond shell crossings 


by using the equation of motion of trajectories, which 
reads as 




dip 
dx ’ 


(A9) 


where p is the rescaled gravitational potential dSD- This 
is the ID version of Ea. (E^ and it explicitly shows that 
particle trajectories obey the same approximate symme¬ 
try, before and after shell crossings. 


2. ID background density perturbation 

To derive the ID consistency relations, we follow the 
method described in the main text for the 3D case, see 
also the Appendix in [2^. As in Ea. d^ . we consider 
two universes with close cosmological parameters, a'{t) = 
a(<)[I — e(<)] and p'{t) = p{t)[l + e(t)]. Substituting into 
the “ID Friedmann equation”, we again find that e(t) 
obeys the same equation as the ID linear growing mode 
D+(t), and we can write e{t) = eoT)+(t). 

Next, the change of frame described in Eg. (1^51) be¬ 
comes 

x' = {1 + e)x, S' = d — e(l -1-5), v' = v + iax, (AlO) 

and at linear order over both 5 and e we have 5l = S'j^ + t. 
This means that the background density perturbation e is 
again absorbed by the change of frame, with eg = 

The redshift-space coordinate s now transforms as 

(All) 

Then, as in Ea. d^ . the redshift-space density contrast in 
the actual unprimed frame, with the uniform overdensity 
A5 lo, writes as 

fc ^ 0 : 5^k,t) = 5*[(1 - e - e/Er)fc,D+,„,/,J,(A12) 

where we disregarded the Dirac factor that does not con¬ 
tribute for wave numbers k ^ Q. Therefore, the derivative 
of the redshift-space density contrast with respect to eg 
reads as 


s' = 1 


d5^{k,t) ^ dD+,, dS^ df,, dS^ 

deo deo dD+ deo df 

-{l + f)D+k^. (A13) 

As shown in the derivative of the linear growing 
mode is dD+/deo = D^, which means that D'^ = D+ + 
egD^. Then, using a' = a — egD^a and the definition 
(IA4I) for / and /', we obtain /' = f + /{e + e/H), whence 

= fD+{l + f). (A14) 

eo=0 

Therefore, Ea. (IA13l) gives Ea. dTTl) . 


dP+eo ^ 
deo deo 
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3. ID consistency relations 


5. Check of the ID consistency relations 


Using the result (ED, the ID version of the consistency 
relations (ED writes as 


±k' 




2=1 


1 + /i 


n d In D+i df. 




n' dk, 

1=1 


(A15) 


The 3D angular average f dUfc//(47r) of Ea. (l37)) is re¬ 
placed by the ID average ^ X]±fc' the two direc¬ 
tions of k' (i.e., the two signs of k'). We again defined 
the reduced polyspectra as in Ea. (l38l) . (d®(fci)..d®(fc„)) = 
{6^{ki)..S^{k„)y SD{kl + .. + kn). 

On large scales we recover the linear theory, with 
S^{k',t') ~ D+(t')(l -|- f')SLo{k'), and Ea. (IA15p also 
writes as 

=PL{k\t') 

^ \ ^ J / k'^0 


n 

^ D' I 

2=1 + ^ 


1 + /i 


n d In D+^ dfi 


-(i+/.)D«5-bi. ® 


1=1 


dkn 


(A16) 


When all times are equal, t' = ti = .. = tn = t, this 
simplifies as Eq. ED- 


4. Zel’dovich solution 

In the ID case, the Zel’dovich approximation is exact 
until shell crossing [2^, . It corresponds to taking for 

the particle trajectories the linear prediction, 


with 


4- 


x{q,t) = q + 'i>L{q,t) 


r+°° Ah 

iq) = ij (A 18 ) 


(A17) 


L[q = 


Therefore, the redshift-space coordinate dMl) writes as 
(using V = ax) 


s = g+(l + /)^L, 


(A19) 


and the redshift-space nonlinear density contrast ED as 

Eq.dZSl). 


a. Impact of a large-scale perturbation on the nonlinear 
redshift-space density contrast 

To check the validity of the ID consistency relations 
from the exact solution (1731) . we simply need the change 
of the nonlinear redshift-space density contrast S^{k) 
when we make a small perturbation Ad^o to the initial 
conditions on much larger scales. Let us consider the 
impact of a small large-scale perturbation A6lo to the 
initial conditions. Here we also restrict to even pertur¬ 
bations, ASLo{—k') = AdLo{k'), as the consistency rela¬ 
tions studied in this paper apply to spherically averaged 
statistics, which correspond to the ±k' averages in the 
ID relations (IA15I) - (IA16|) . Then, expanding Eg. (1751) up 
to first order over ASlq, and over powers of fc', we obtain 


k'^0: AS%k) = {1 + f)D+ 


dk'Ahoik') 


X / ^ e-''k,+kii+f)D^ f ^hoik'd ^ikq).iA20) 

J Stt 

Here the limit fc' —>■ 0 means that we consider a pertur¬ 
bation of the initial conditions AdLo{k') that is restricted 
to low wave numbers, k' < A, with a cutoff A that goes 
to zero (i.e., that is much smaller than the wave numbers 
k and 27r/g of interest). 

On the other hand, from the expression (l73ll we obtain 
at once the exact result 


27T 


kq 


dS^ d6^ 5 ( 5 ® 

X ^k{l+f)D+ f f)(ikq). (A21) 

The comparison with Eq. (IA20I) gives 

' dS^(k) 


k'-rO: A5%k)=D+ 


dk'ASLoik’) 


dhiD 


-I- 


+ + .(A22) 


01 


dk 


The consistency relations (IA15p - (IA16l) and (172]) only 
rely on the expression ED, which also reads (at linear 
order over cq) as 


51nD 


+ 


df 


-(l + /)fc 


dS^jk) 

dk 


(A23) 


Since we have eg = ASlo = f dk'ASi,o(k'), we recover 
Ea. (IA22p . This provides an explicit check of Eq. ED, 
hence of the ID consistency relations. 
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b. Explicit check on the redshift-space density polyspectra 

Instead of looking for the impact of a large-scale lin¬ 
ear perturbation on the nonlinear density contrast, as in 
Sec. lAS^ we can directly check the consistency relations 
in their forms (IA15I) or (1721) . Considering for simplic¬ 
ity the equal-time polyspectra (17^ . we define the mixed 
polyspectra, formed by one linear density contrast and n 
nonlinear redshift-space density contrasts. 


.. -|- kn-i)- Finally, in the limit fc' —>■ 0 we expand the 
terms up to first order over fc', and we obtain 

Ei = PUm + f) J 

g-i kj qj 

X g-(l-H/)='£>+/2/dfe/fc" 

Proceeding in the same fashion, the n—point redshift- 
space polyspectra read as 


n—1 

1 + i X! 


i=i 


X g(l+/)-D+ /dfc/fc &Lo(.k) E”=i fei 

where in the last expression we used Eo. dTSl) . The Gaus¬ 
sian average over the initial conditions Slq gives 



P:={S^ik,,t)J%kr.,t)y 

= f dgi..dg„-i ijy:-fk,q, 

J (27r)"-i 

X g-(i-H/)"Di/2/dfe/fc= Pi.o(fc)|E"Ti'feAe‘'''’^-i)r,(A27) 


El 


PUk') 

k' 


(1 + /) 


dgi..dg„ 

(27r)’" 






(A25) 


Making the changes of variable qi = q'^ qn, ■■■, q-n-i = 
q'n-i+qn, the argument of the last exponential does not 
depend on qn- Then, the integration over qn yields a 
Dirac factor Soik' ki kn), that we factor out by 

defining E^ = E^Snik' -|- fci -I- .. -f A:„), with a primed 
notation as in Ea. (l551) . and we replace kn by —{k' -\-ki-\- 


Then, we can explicitly check from the comparison with 
Ea. (IA26p that we have the relation 


k'^O: Ei = PL{k') 




-(1 + /) 


n —1 


d 


ps 


(A28) 


d In ki 

i—1 

and we recover the consistency relation GH). [In 
Eq. (IA28P the right-hand side does not involve because 
it has been replaced by —{ki + .. + fc„_i) in Ea. (IA27p . 
using the Dirac factor Soiki -I- .. -I- kn).] 
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